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ABSTRACT 

We study the phase space available to the local stellar distribution using a Galactic potential con- 
sistent with several recent observational constraints. We find that the induced phase space structure 
has several observable consequences. The spiral arm contribution to the kinematic structure in the 
solar neighborhood may be as important as the one produced by the Galactic bar. We suggest that 
some of the stellar kinematic groups in the solar neighborhood, like the Hercules structure and the 
kinematic branches, can be created by the dynamical resonances of self-gravitating spiral arms and 
not exclusively by the Galactic bar. A structure coincident with the Arcturus kinematic group is 
developed when a hot stellar disk population is considered, which introduces a new perspective on the 
interpretation of its extragalactic origin. A bar-related resonant mechanism can modify this kinematic 
structure. We show that particles in the dark matter disk-like structure predicted by recent LCDM 
galaxy formation experiments, with similar kinematics to the thick disk, are affected by the same 
resonances, developing phase space structures or dark kinematic groups that are independent of the 
Galaxy assembly history and substructure abundance. We discuss the possibility of using the stellar 
phase space groups as constraints to non-axisymmetric models of the Milky Way structure. 
Subject headings: dark matter — Galaxy: disk — Galaxy: evolution — Galaxy: kinematics and 
dynamics — Galaxy: structure — solar neighborhood 



1. INTRODUCTION 

The study of stellar kinematic groups in the solar 
neighborhood has a long tradition in Galactic astron- 
omy, going ba ck to the disco very of the Hyades and U rsa 
Major groups l|Proctoilll86"9h . In lAntoja et~al] (|2008f ) we 
presented a study of the solar neighborhood kinematic 
groups using a sample of 24,190 stars. We confirmed the 
existence of the Sirius, Coma Berenices, Hyades-Pleiades 
and Hercules branches (Fig. [Hi). They all have a negative 
slope of 16° in the U-V plane (following the standard 
definitions of U and V) and a slight curvature. 

Considerable amount of work has been performed in 
an attempt to explain the origin of these kinematic 
groups. External pr ocesses like past accretion events 
(|Navarro et al.l I2004D and internal disk processes like 
star formation bursts or secular dynamics have been 
proposed. Although they were initially considered mu- 
tually exclusive, all such mechanisms are natural in 
current galaxy formation models ([ Klypin ct al. 2008; 
lRomano-I)i'az et al.l 120081: iCeverino fc Klvpin ,2007b ). 
Recently, the hypothesis of the disk-dynamical ori- 
gin of these structures has gained popularity par- 
tially because of the consistency of the Hercules struc- 
ture with the effects o f the Galac t ic ba r resonances 
(iKalnaid fl99l iDehnenI [2000t [200l . In addi- 

tion, steady or transient spiral arms were proposed 
to explain the charact eristics of s ome of t hese stel- 
lar groups (iSkuli a n et al. 1999: De Sirnone et al. 2003; 
iFamaev et al.l 120051: lOuillen fc MinchevI l2005r . A more 
recent study argues that the combined effect of a bar and 
spiral arms is necessary to accur ately reproduce obse rva- 
tions in the solar neig hborhood (|Chakrabartvll2007D . In 
that case, spiral arms at the weaker end of the observa- 
tional range (with a fractional amplitude of less than 4% 



of the background disk density) were used. However, in 
that study the spiral arms seem to contribute only in the 
fine structure, thus weakening constraints on the Milky 
Way (MW) spiral structure based on the solar neighbor- 
hood kinematics. Previous studies have considered mod- 
els for the non-axisymmetric components of the Galaxy 
motivated by the dynamical point of view, in particu- 
lar weak spiral arms described by a cosine function and 
bars described by the quadrupole perturbation. How- 
ever, it is unclear whether there is any dependence of 
the induced local solar neighborhood kinematics on the 
detailed Galactic structure. Moreover, the initial condi- 
tions hardly consider the evolution of the MW. 

Our contribution improves upon earlier studies. First 
we investigate the stellar kinematic response to a model 
that satisfies published observational constraints to the 
MW structure, in particular to the non-axisymmetric 
components. Secondly, we consider for the first time 
initial conditions and integration times that attempt to 
represent stars born at different times and with different 
kinematic conditions, like those in the solar neighbor- 
hood. Lastly, we investigate effects on the local dark 
matter kinematics, in particular in the disk-like dark 
matter structure recently predicted by LCDM models. 
This issue is important to predict signals in direct dark 
matter detection experiments. 

2. SIMULATIONS 

In order to study the effect of the non-axisymmetric 
Galactic structure on the solar neighborhood kinematic 
distribution, we have performed numerical integrations 
of test particle orbits on the Galactic plane, adopting 
the initial conditions discussed in Sect. 12.21 and the po- 
tential described in Sect. 12.11 Each particle integration 
time is initialized at a value t = — r and ended at t= 0, 
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TABLE 1 

NON-AXIS YMMETRIC GALACTIC DiSK CO MPONENTS 
(jPlCHARDO ET AL.|[200l [200ll 



Parameter 


Value 


Bar axis ratio 


10:3.12 


Bar scale lengths ( kpc) 


1.7 and 0.54 


Bar angle respect to the Sun (°) 


20 


Bar mass (Mq) 


lOlO 


Bar pattern speed (see text) fij, (kms~^kpc~^) 


60, 45 


Spiral arm pitch angle (°) 


15.5 


Spiral arm scale-length ( kpc) 


2.5 


Spiral arm mean force ratio (%) 


10 


Spiral arm pattern speed f2sp (kms^^kpc"^) 


20 



where r is the time during which the particle is exposed 
to the non-axisymmetric potential. The induced kine- 
matic distribution at the end of the simulation is studied 
by considering the particles inside a circle of radius 500 pc 
centered at the solar position. Finally, the predicted and 
the observed distributions are compared. The motion 
equation s were integrated with the Bulirsch-Stoer algo- 
rithm of lPress et al.l ()2004l ). conserving Jacobi's integral 
within a relative variation of \{Eji — Ejf)/ Eji\ k, 10"^^ 
for only-arms and only-bar models. The reference frame 
used for the calculations is the rotation frame of the spi- 
ral arms when only this non-axisymmetric component is 
considered and the one of the bar in the other cases. In 
all cases, we check that the number of particles in the 
final distributions is statistically robust. 

2.1. The Galactic Model 

We use th e Ga lactic potential described by 
iPichardo et al.l ()2003l . 2004) consistent with several 
recent observational constraints. For our study the most 
relevant parts are the bar and spiral arm components, 
whose corresponding parameters are presented in Table 
[TJ The model is compared with observation s and 
with other models elsewhere (jPichardo et al.l l2003l . 
|2004[ ). The bar consists of a prolate mass distri bution 
that resembles Model S of iFreudenreichI (|1998f ) from 
COBE/DIRBE data. We consider two alternative bar 
pattern speeds (60 an d 45 kms~^kpc ~ i^) ac cording 
to the range found by iDebattista et al.l ()2002[ ). The 
spiral arm model is a 3D steady m odel with two arms 
that traces the locus reported by iDrimmel &: Spergel 
(|200ll) using K-band observations, and also in agreement 
with the latest results from Spitzer (Benjamin 2009, 
private communication). The solar radius is close to the 
bar outer Lindblad resonance (OLR) and to 4:1 inner 
Lindblad resonance (ILR) of the spiral arms. The ratio 
of the radial force of the spiral arms to the background 
peaks at '^1.5:10 and the mea n value along t he arm is 
approximately 1:10 (see Picha rdo et al.|[2003f). This is 
at the higher end of the hmits suggested bv lPatsis et al.l 
(|1991h for MW type galaxies. Instead of introducing 
the spiral arms and bar adiabatically to avoid transient 
features on the model, we prefer to test self-consistency 
through an analysis o f the stellar orbital reinforcement 
of the potentia l , as in iPatsis et al.l ()199lD . presented in 
IPichardo et all ((2001 [200ll 

2.2. Initial Conditions 

The choice of appropriate initial conditions is some- 
what controversial. Traditionally they are motivated by 



self-consistency with the present stage of the Galactic 
disk structure. However, current scenarios of galaxy for- 
mation and evolution predict that both external and in- 
ternal perturbation mechanisms affect the disk kinemat- 
ics. The integration time in our simulations is intended 
to consider only the latest stages in the Galaxy evolu- 
tion. This late evolution was aff ected or even domi nated 
by secular dynamical processes ()Klvpin et al.ll2008l ). We 
avoided longer integration times (more than 2 Gyr) be- 
cause they might require the inclusion of Galactic evo- 
lution like bar weakening or transient arms, as well as 
the effect of external perturbations to the MW disk, 
thus masking the effect of the current MW structure. 
Therefore we focused on the recently induced kinematic 
structure in the solar neighborhood. We explored three 
different types of initial conditions (ICl, IC2, ICS), all 
of them assuming an exponential disk scale length with 
Rp = 2.5 kpc. 

ICl — The initial velocity distribution relative to the 
Regional Standard of Rest is adopted as a Gaussian (with 
(7(7 = uv = 5 km s~^) constant for all radii. Each particle 
is exposed to the non-axisymmetric perturbations for a 
time (r) chosen at random between and 2 Gyr. This 
maximum integration time corresponds approximately to 
14-20 and 6 revolutions of the bar and the spiral arms, 
respectively. With this set of initial conditions we aim to 
simulate the birth of stars in the disk with small velocity 
dispersion and the effect of having stars with different 
ages (different integration time) at the solar position^ . 

IC2 — The distribution function satisfies the coUision- 
less B oltzmann equation as it is discussed in Hcrnquisl 
(|1993f ). The velocity dispersion has an exponential pro- 
file with scale length R„ = 7.5 kpc and local normal- 
ization tT(i?0 = 8.5 kpc) ^ 20kms~^. In this case the 
integration time is fixed at r = 400 Myr for all particles 
(corresponding approx imately t o 3-4 bar revolutions and 
similar to the value of iDehnenl ([2000)). With these ini- 
tial conditions, we can study the relatively rapid induced 
effects of the non-axisymmetric component on the local 
kinematics. 

ICS — Identical to IC2 but with a higher velocity dis- 
persion of cr(i?0 — 8.5 kpc) ^ 40kms~^, closer to an old 
thin disk component, or to the thick disk. Recent stud- 
ies of the formation of MW mass galaxies predict that a 
flattened dark matter structure mirroring t he thick disk 
properties will form in a LCDM Universe (jBruch et al.l 
120081 : iRead et al.ll2008l ). Therefore, these initial condi- 
tions are also consistent with particles in the dark disk. 

3. SPIRAL ARM AND BAR CONTRIBUTION TO THE 
LOCAL KINEMATIC STRUCTURE 

Models Using Only Spiral Arms — In this section we 
present the kinematic structure developed by the self- 
gravitating spiral arm model. We find that this model 
produces ample substructure. In particular, assuming 
low velocity dispersions (ICl), it reproduces a branch at 
low angular momentum (with V ^ — 40kms~^ consis- 
tent with the Hercules structure) separated from a cen- 
tral group of substructures. To our knowledge, it is the 

^ N-body disk simu lations assuming a low Q parameter (e.g. 
ITho masson et al."199in show that a configuration with strong spi- 
ral arms and low velocity dispersion can be sustained for at least 7 
disk revolutions. These experiments are similar to our integrations 
using ICl. 
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Fig. 1. — Density field in the U-V plane obtained by wavelet denoising for: a) the observed sample; b, c) simulations with a model using 
only spiral arms; d, e, f, g) simulations with a model using only bar; h, i) simulations with a combined model (bar + spiral arms). Notice 
the difl'erent scales of the panels. 

first time that an unbarred model has produced a sim- 
ilar structure. Figure [TJa shows the U-V plane for the 
solar position. Although the shape of the structure at 
V = — 40kms~^ is not exactly equal to the Hercules 
branch, it does show that the spiral arms by themselves 
crowd the velocity space at these negative V. In par- 
ticular, the test particles in this kinematic group have 
a perturbation exposure time of at least 1 Gyr. The 
central part of the distribution seems to be split into 
two groups or branches, resembling some of the observed 
central kinematic branches (Fig. [1^). In order to ex- 
plore small variations re spect to the configuration of 
iDrimmel fc Spergell (|2001[ ). given the uncertainty on the 
exact position of the arms, we also show in Figure [!}: 
the U-V plane for a region located at </> = —40°, that is 
at the solar circle but at 40° of the Sun in the counter- 
clockwise direction. In this case and for a wide range of 
0, the equivalent to the Hercules structure is also gen- 
erated. We conclude that the contribution of the spiral 
arms to the solar neighborhood kinematics may be com- 
parable to that of the bar. The sensitivity of our results 



to the properties of the arms indicates that local kinemat- 
ics can be used as one of the constraints to the current 
observational ambiguities about this non-axisymmetric 
component. 

Models Usi ng Only a Bar — In agreement with previous 
studies (e.g. lDehnenll2000f ). we find that a Galactic model 
using only a bar can trigger a kinematic group at neg- 
ative V and U resembling the Hercules branch. How- 
ever, in our experiments this branch can be populated 
only under some combinations of model parameters and 
initial conditions. For example, a bar pattern speed of 
= eOkms-^kpc^^ [Rolr/R = 0.74) under IC2 pro- 
duces a branch at an incorrect position ([/ > 0). How- 
ever, if the imaginary solar observer is located at an inner 
radius i? = 6 — 7kpc (where Rqlr/R = 0.97) a struc- 
ture appears at V ^ — 35kms~^ and [/ < 0, as is shown 
at the top left of Figure [TJi. A kinematic group at V, 
U < also appears if the bar pattern speed is reduced to 
45kms^^kpc~^ (Fig-Hfe), although the mean branch ve- 
locity V is now slightly smaller: V ^ — 20kms^^ (again 
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with Rolr/R = 0.97). In summary, the Hercules branch 
at the correct position may be achieved by varying either 
the model pattern speed, the observer position or even 
the initial conditions. 

If we consider hotter initial conditions (ICS), arch- 
shaped structures appear at lower V. These struc- 
tures appear even if we consider an axisymmetr i c disk 
model, as also recently reported bv lMinchev et al.l ()2009D 
who assumed a non-relaxed model. In our simula- 
tions the positions of these kinematic arches are mod- 
ified when the bar is added to the model. This sug- 
gests that both kinematic initial conditions and Galac- 
tic structure contribute to create the kinematic groups. 
For example, differences are observed if we use either 
Qb = 45kms-ikpc"^ (Fig. [If ) or ilt = eOkms-^kpc"^ 
(Fig.[T^). The arches developed a,t V ^ — 40kms^^ and 
— lOOkms"^ for the lower bar pattern speed (not seen 
when Qij = 60kms~^kpc~ ) are closer to the Hercules 
and Arcturus observed structures. ^ This result should 
be considered only as evidence of sensitivity to the model 
but not necessarily as favoring a particular Q,b value. Fur- 
thermore, these simulations show the important role of 
the bar in the development of the local kinematic struc- 
ture. 

Combined Models — Here we present the results obtained 
with the model that includes both the bar and the spi- 
ral arm perturbations. This model, with different pat- 
tern speeds for the bar and the spiral s, is supported by 
recent studies like iPatsis et al.l ()2009t ). They reported 
that the best match between observations and models 
for the galaxy NGC SS59 is found when different pattern 
speed for the bar (39 kms^^kpc^^) and for the spiral 
arms (15kms^^kpc^ ) is considered. Now we show ex- 
amples of how the branches generated by the spiral arms 
or the bar are maintained when the combined model is 
used. Figure[T)i shows the results for the spiral arms plus 
the bar with fif, = 45kms~^kpc~^ used with ICl near 
the solar position (0 — —40°). In Figure [TJ, we show 
the case corresponding to spiral arms and a bar with 
fib = 45kms~^kpc^^ under ICS. Here the arch-shaped 
structures associated with Hercules and Arcturus are still 
observed. 

The diversity of the initial conditions in our simula- 
tions is not extensive. Moreover, the observed velocity 
field is likely to be the result of the combined action of 
several Galactic processes. Therefore, comparison be- 
tween the observed velocity distribution^ (Fig. [I^) and 
the results of test particle simulations is not straightfor- 
ward. A first approach might be to compare the observed 
velocity field with the combined model results under ICl 
and ICS, where the central branches and the kinematic 
groups at low angular momentum would appear simul- 
taneously. The next complexity level might require the 
inclusion of age and kinematic criteria in the definition of 
the artificial observational sample. We defer such details 
to a future study, since they do not alter the conclusions 
of this Letter. 

4. LOCAL DARK MATTER KINEMATICS 

^ Notice that the observed Arcturus structure is out of the range 
of Figure [iti as t his structure is not promi nent in the sample u sed 
bv lAntoiaerall ll2008i '). See Figure 2 of iWilliams et al.l (120081) for 
a U—V plane containing the Arcturus structure. 



Another unexpected aspect of the bar- and spiral arm- 
induced phase space structure is the effect on the lo- 
cal dark mat ter kinematics. P articles in a possible dark 
matt e r bar (IColm et al.l [20061 ) . dark disk (jBruch et al.l 
120081 [Read et al.' '2008") and even in the dark matter 
halo fAthanassoula 2005) are trapped/scattered in the 
same resonances as stars are. It has been shown that 
different bar-induced resonai ices can be popula, t ed by 
disk-like and halo parti cles (jAthanassoulal l2002l I200SI : 
iCeverino fc Klvpinll200l . 

Our stu dy considers only the case of the dark disk. 
Following l|Read et al.ll2009| ) we assume that our ICS ini- 
tial conditions set also represents the dark particles in 
the dark thick disk. Thus Figures [Tf, [1^ and[TJ, obtained 
from ICS, also reproduce the local dark matter kinemat- 
ics induced by our MW dynamical models. Our results 
show that these models generate dark matter currents in- 
side the Galactic dark disk. These dark-matter currents 
would be independent of the Galactic assembly history 
or the dark substructure abundance. The dynamical his- 
tory of the Galaxy and its detailed large scale structure 
may help to establish whether the amplitude of the dark 
matter kinematic structure is detectable by planned dark 
matter detection experiments. 

5. DISCUSSION AND CONCLUSIONS 

The orbital analysis of a MW model consistent with 
published observational constraints has several unex- 
pected consequences for the local stellar kinematics. 

The spiral arm contribution to the resonant structure 
in the solar neighborhood may be comparable to that 
of the Galactic bar. The main differences to previous 
studies are the arm force contrast and force field shape 
( Pichardo et al.. 200S. Fig. 5, 8 and 9), as well as the va- 
riety in initial conditions. In particular, we find that the 
Hercules structure may be produced by the spiral arms 
and not e xclusivel y by b ar resonances as traditionally 
believed. iD ehnenI (|2000( ) concluded that the Hercules 
branch is unlikely to have been produced by resonant 
scattering processes due to spiral arms. His main argu- 
ment is that spiral arms do not act on stars with epicy- 
cles greater than the interarm separation. However, the 
epicycle amplitude of the stars in the structure created in 
our simulations^ (at U ~ [-30, SO] and V ~ [-45, -35]) 
is about Skpc. On the other hand, the interarm sepa- 
ration in the different loci of our model is in the range 
5.5 — 7kpc. This is sufficient to prod uce the Hercule s 
structure and still in agreement with Dchnejil (|2000( ). 
Other authors like Quillen & Minchev (2005) reproduced 
the Hyades-Pleiades and Coma Berenices branches as as- 
cribed to periodic orbits related to the spiral arm 4:1 ILR, 
but not the Hercules structure. For some parameter com- 
binations our model reproduces both kinds of structure. 
Detailed comparison with other authors is not straight- 
forward because many of them used the tight-winding 
approximation to model the spiral arms, a different sim- 
ulation strateg}^j_MdJn some cases, even a four-armed 
model (|Chakral3artvll2007[ ). while ours includes two arms. 

As in previous studies, we confirm that the Hercules 
structure can also be produced by the Galactic bar under 
certain combinations of the model parameters and initial 

^ This val ue is obtained the oretically from the epicycle approxi- 
mation (e.g. lAsiain et al.ll999l) or directly from orbital integration. 
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conditions. However, we considere d a prolate bar (for the 
detailed force field see Fig. 16 of lPichardo et al.ll2004D 
instead of the widely used cosine bar potential (|Dehnenl 
120001) . The study of possible differences in the local kine- 
matic structure generated by each model is postponed for 
a future study. 

A subset of our experiments develops a structure re- 
sembling the Arcturus kinematic group. The required 
condition seems to be a relatively hot stellar disk popu- 
lation similar to the thick disk in kinematics. This kine- 
matic grou p may have arisen from a past accret ion event 
(jNavarro e t al. 2004; Villal obos fc Helmil l2009f ) . Alter- 
natively IWilliams et al.. (,2008h recently discussed a pos- 
sible disk-dynamical origin of Arcturus based on stel- 
lar population evidence, and postulated the bar 6:1 
OLR as the tr i ggerin g mechanism. On the other hand, 
iMinchev et al.l ()2009D proposed an origin related to non- 
equilibrium initial conditions. Our results support an 
internal disk origin and, in addition, we find that the 
dynamics of the bar has an influence on these low an- 
gular momentum kinematic groups as it modifies these 
structures generated in the axisymmetric model. Pre- 
liminary results of orbit integration in 3-D show that the 
main arches generated in the 2-D case are maintained and 
the Galactic bar significantly takes part in defining their 
shape. However, the 3-D case deserves special attention 
and is the subject of ongoing studies. The fact that a 
vrms of 40 km s~^ is high enough to allow the creation of 
these arches in our experiments is encouraging, but it is 
still not sufficient to disentangle the origin of Arcturus. 
A deeper discussion of this point is beyond of the scope 
of this letter, but it is the subject of an upcoming study. 

We show that the Galactic non-axisymmetric poten- 
tial develops dark kinematic groups in the dark disk pre- 



dicted in cosmological simulations of galaxy formation. 
These currents are independent of the halo substructure 
abundance and may have new consequences for planned 
dark matter direct detection experiments. 

The dependence of the stellar kinematics in the so- 
lar neighborhood on the structure, dynamics and ini- 
tial conditions of our experiments suggests that kine- 
matic groups may provide a useful constraint on non- 
axisymmetric MW models. A systematic scan of the 
associated parameter space will be required in order to 
disentangle the origin of the different kinematic groups 
in the solar neighborhood. We are currently exploring 
strategies to do so. Studies like the ones discussed in 
this letter will derive benefit from upcoming surveys like 
GAIA and SEGUE2. In summary, the imprints of the 
non-axisymmetric Galactic structure on the local stellar 
kinematics are strong. 
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